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Abstract. Magnetic fields correlated on kiloparsec scales are seen in spiral 
galaxies. Their origin could be due to amplification of a small seed field by 
a turbulent galactic dynamo. We review the current status of the galactic 
dynamo, especially the constraints imposed by magnetic helicity conservation. 
We estimate the minimal strength of the large-scale magnetic field which could 
arise inspite of the helicity constraint. 
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1. The galactic dynamo 

Magnetic fields in spiral galaxies have strengths of order few lO^^G, and are coherent 
on scales of several kpc (Beck et al. 1996). In several disk galaxies, like M51 and NGC 
6946, they are also highly correlated (or anti-correlated) with the optical spiral arms. 
How do such ordered, large-scale fields arise? One possibility is the dynamo amplification 
of a weak but nonzero seed field. We critically review here the operation of the galactic 
dynamo, particularly emphasing the constraints which arise due to the conservation of 
magnetic helicity in highly conducting plasma. 

The evolution of the magnetic field is described by the induction equation 

— = V X (v X B -77V X B). (1) 
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Here B is the magnetic field, v the velocity of the fluid and ry the resistivity. B = is 
a perfectly valid solution of the induction equation. So there would be no magnetic field 
generated if one were to start with a zero magnetic field. There arc a number of battery 
mechanisms, invoking small additional source terms to Ohms law, which lead to a seed 
magnetic field from zero fields (cf. Rees 1994; Subramanian, Narasimha and Chitre 1994). 
This seed field is gcncrically much smaller than the galactic fields. Therefore some form 
of dynamo action, due to motions which act to exponentiate small seed fields efficiently, 
is essential to explain observed galactic fields. 

Galactic dynamos depend on the following two features: First, disk galaxies are 
differentially rotating systems. Also the magnetic flux is to a large extent frozen into the 
fluid. So any radial component of the magnetic field will be efficiently wound up and 
amplified to produce a toroidal component. But this results in only a linear amplification 
of the field. To obtain the observed galactic fields starting from small seed fields one 
needs a way to generate the radial component from the toroidal one. If this can be done, 
the field can grow exponentially and one has a dynamo. 

A mechanism to produce the radial field from the toroidal field was invented by Parker 
(1955), and is known as the a-effect (Steenbeck, Krause and Radler 1966). The essential 
feature is to invoke the effects of cyclonic turbulence in the galactic gas (cf. Ferriere 
1998). The interstellar medium is assunicid to be turbulent, due to for example the effect 
of supernovae randomly going off in different regions. In a rotating, stratified (in density 
and pressure) medium like a disk galaxy, such turbulence becomes helical. An upward 
moving fiuid parcel, expands and the coriolis force makes it rotate retrogade, generating 
negative kinetic helicity in the northern hemisphere. Downward moving fluid contracts, 
and the coriolis force now makes it rotate in the prograde direction. This contributes to 
helicity of the same sense. Helical motions of the gas perpendicular to the disk draws 
out the toroidal field into a loop which looks like a twisted f2. Such a twisted loop is 
connected to a current which has a component parallel to the original toroidal field. If 
the motions have a non-zero net helicity, this parallel component of the current adds 
up coherently. A toroidal current then results from the toroidal field. Hence, poloidal 
fields can be generated from toroidal ones. (Of course microscopic diffusion is essential to 
make permanent changes in the field). This closes the toroidal-poloidal cycle and leads 
to exponential growth of the mean field. 

In quantitative terms, suppose the velocity field is the sum of a mean, large-scale 
velocity Vq and a turbulent, stochastic velocity vy. The induction equation becomes 
a stochastic partial differential equation. Split the magnetic field B = Bq -I- b, into a 
mean field Bq =< B > and a fiuctuating component b. Here the average <>, is defined 
either as a spatial average over scales larger than the turbulent eddy scales (but smaller 
than the system size) or as an ensemble average. Assume the turbulence to be isotropic, 
homogeneous, helical and have a short (ideally delta function) correlation time r. Then 
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one can derive the mean- field dynamo equation for Bq, 

= V X (Vo X Bo + £ - r?V X Bo) . (2) 

£ =< vt X b >= aoBo - Vt^ x Bo- (3) 

Here e is the turbulent cmf, = — (t/S) < V7-.(V x v^) >, is the dynamo a-effect, 
proportional to the kinetic helicity and rjT ~ {t/ 3) < v"^ >, is the turbulent magnetic 
diffusivity proportional to the kinetic energy of the turbulence. This kinematic mean-field 
dynamo equation, has exponentially growing solutions, provided a dimensionless dynamo 
number has magnitude D = \aoGh^ri^'^\ > Dcrit ~ 6 (Ruzmaikin, Shukurov and SokolofF 
1988). (Here h is the disk scale height and G the galactic shear, and we have defined D 
to be positive). While the a-effect is crucial for regeneration of poloidal from toroidal 
fields, the turbulent diffusion turns out to be also essential for allowing changes in the 
mean field fiux. The mean field grows typically on time-scales a few times the rotation 
time scales, of order 10^ yr. Modulations of a, and rjT due to the spiral arms, can also 
lead to large-scale fields, correlated (or anti-correlated) with the optical spirals (Mestel 
and Subramanian 1991; Moss 1998; Shukurov 1998). 

The kinematic mean-field equation neglects the the back-reaction on the velocity due 
to the Lorentz forces. This rapidly becomes a bad approximation, due to the more rapid 
build up of magnetic noise compared to the mean field (Kulsrud and Anderson 1992). 
Both direct numerical simulations of the non-linear dynamo (Brandenburg 2001, Branden- 
burg and Sarson 2002, Brandenburg, Dobler and Subramanian 2002) and semi-analytic 
modelling of the non-linear effects (Subramanian 1999; Brandenburg and Subramanian 
2000) point to the crucial role played by magnetic helicity conservation in limiting mean 
field growth. 



2. Magnetic helicity conservation and the galactic dynamo 

The magnetic helcity associated with a field B = V x A is defined as H = J A.B dV, 
where A is the vector potential (Moffat 1978, Berger and Field 1984). Note that this 
definition of helicity is only gauge invaraiant (and hence meaningful) if the domain of 
integration is periodic, infinite or has a boundary where the normal component of the 
field vanishes. In this case, under a gauge transformation A ^ A — V^", the additional 
term in the helicity, / V?/'.B = / tp'B.dS — / tpV.'B dV = 0. H measures the linkages 
and twists in the magnetic field. From the induction equation one can easily derive the 
helicity conservation equation, 

f = -2,/^J.B.F, (4) 

where J = (c/47r)V x B is the current density. So in ideal MHD with 77 = 0, magnetic 
helicity is strictly conserved. However, this does not guarantee conservation of H in the 
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limit — *■ 0, because the current helicity, / J.B dV , may still become large. For example, 
the Ohmic dissipation rate of magnetic energy Qjouie = r]{^'K/(?) J J'^dV can be finite 
and balance magnetic energy input by motions, even when r) ^ 0. This is because small 
enough scales develop in the field (current sheets) where the current density increases 
with decreasing rj as (x ry"^/^ as 77 — > 0, whilst the rms magnetic field strength, firms, 
remains essentially independent of rj. Even in this case, however, the rate of magnetic 
helicity dissipation decreases with rj like oc 77"*"^/^ — > 0, as ?7 ^ 0. Thus, under many 
astrophysical conditions where Rm is large {t] small), the magnetic helicity H, is almost 
independent of time, even when the magnetic energy is dissipated at finite rates. 

Coming back to the mean-field dynamo, we note that its operation automatically leads 

to the growth of linkages between the toroidal and poloidal mean fields. Such linkages 
measure the helicity associated with the mean field. One then wonders how this mean field 
(galactic) helicity arises? To understand this, we need to split the helicity conservation 
equation into evolution equations of the sub-helicities associated with the mean field, say 
Hq = J Aq.Bq dV and the fluctuating field h = J < a.b > dV =< a.b > V. The 
evolution equations for Hq and h are 



Here, and henceforth, we assume that the surface terms can either be neglected or are 
zero (becuase of boundry conditions). We see that the turbulent emf e transfers helicity 
between large and small scales; it puts equal and opposite amounts of helicity into the 
mean field and the small-scale field, conserving the total helicity H = Hq + h. So if one 
were to start with zero total helicity intially, in a system with large Rm, one will always 
have i?o + /i ~ 0, or \Ho\ w \h\. 

Note that for a given amount of helicity, the energy associated with the field is in- 
versely proportional to the scale over which the field varies. If for example, the small-scale 
field were maximally helical, and varied on a single scale, with associated wave number, 
kf, we will have kf < a.b >=< b^ >. Similarly in a periodic box, a maximally helical 
large scale field with wave number km, satisfies, km J Aq.Bq = / dVBQ. (Hence- 
forth, we will denote the volume average of mean field quantities Xq over the scale of 
the system, J{dV/V)Xo, by Xq ). So, helicity conservation, with \Hq\ x \h\, implies 
Bq « (km/kf) < b^ >. Now in general, < b^ >^/^, will saturate near the equipar- 
tition field strength, say B^^ = Airpv^. (Here p is the fluid density). So one obtains 
Bg « {km/kf)Bf,g <§; Bgq, for km/kf -C 1. The mean field is expected to attain atmost 
sub-equipartition values for Rm S> 1, if helicity is strictly conserved. 







The galactic dynamo also involves shear and the generation of the toroidal field by 
shear, does not involve the generation of net helicity. A periodic box simulation with an 
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imposed periodic shear (Brandenburg, Bigazzi and Subramanian 2001), suggests that the 
hehcity constraint apphes now to the product of the mean toroidal {Bt ) and poloidal fields 
(Bp). So BtBp/kjn ~ I < a.b > \ ^ S < > /kf. Here 6 < 1 takes into account that the 
small-scale field will also not be fully helical. In galaxies one usually has Bt/Bp = Q > 1. 
This implies Bf ^ {Q5){km/kf)Bl^. In principal one can have large Bt at the cost of 
Bp, with strong shear. These estimates give upper limits to the mean- field strength with 
and without shear. However, in both cases, the limits are smaller than Beq only by the 
square root of the ratio of small to large scales (and by a further factor ((3(5)^/^, in case of 
shear). Whether such mean field strengths are indeed realised, depends also on detailed 
dynamics of mean-field dynamo saturation, to which we now turn. 



Modelling dynamo saturation and minimal mean galactic 

fields 



As a crude model of how the dynamo saturates, when the dynamo is not too supercritical 
(see below), one may use the quasi-linear theory applicable to weak mean fields (cf. 
Pouquet, Frisch and Leorat 1976; Zeldovich, Ruzmaikin and Sokoloff 1983; Gruzinov 
and Diamond 1994; Bhattacharjee and Yuan 1995; Subramanian 2002a). This gives a 
re-normalised turbulent emf, with a = -\- aM, where um — {t/Z) < b.V x b > 
/{AiTp), is proportional to the small-scale current helicity. The turbulent difusion rjT is 
left unchanged to the lowest order, although to the next order there arises a non-linear 
hyperdiffusive correction to e (Subramanian 2002a). One can now look for a combined 
steady state solution to the helicity conservation equation (^) , and the mean-field dynamo 
equation. This work is in progress (Subramanian 2000b), and preliminary results are 
reported here. (Similar work is also being done by Brandenburg and Blackman (private 
communication); see also Brandenburg 2002; Field and Blackman 2002). Assume again 
that the small-scale field has a scale kj^. We then have < b.V x b >= k^ < a.b > 
(although, for fields which are not maximally helical, these helicities can not be related 
to the energy). Using this, we can write dh/dt in Eq. (^, in terms of daAi/dt and hence 
in terms of da/dt. Further, for Rm large enough that helicity is conserved, and for a 
large-scale field which varies on scale fc~^, we can approximately write, (V x Bo).Bo ~ 
fc^Ao.Bo ~ — fc^ < a.b >= —{km/kfY < b.V x b >. This relation is of course strictly 
valid only if Hq is well defined. This inturn requires that negligible mean magnetic flux 
leave the disk, which is likely in thin disk dynamos (cf. Ruzmaikin, Shukurov and Sokoloff 
1988). (A more careful treatment will involve using the relative helicity of Berger and 
Field 1984). The helicity conservation equation then gives a dynamical equation for 
a- quenching 

rirkj dt Bj^ kj r/r 

We see that non-linear effects leads to a decrease in a with time, till the RHS side Eq. (0) 
becomes zero. (Such quenching of a has earlier been discussed by Klceorin and Ruzmaikin 
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1982; Zeldovich, Ruzmaikin and SokolofF 1983). As a decreases, the effective dynamo 
number of the galactic dynamo, D = jaG/i^r/^^l, will also decrease from an initial value 
Do — |aoG/i^?7^^| and lead to a saturation of the mean field growth when D = Dcrit- This 
happens when a = agat = Oif){Dcrit/ Dq). The stationary solution for both the dynamical 
a quenching equation and the mean field dynamo equations is then obtained by equating 
the RHS of Eq. (@) to zero, and substituting a value of a = Usat given above. For 
"HMt <C (km/kf)^, this gives an estimated mean field strength Bmean = |Bo| 

km 

Bmean ~ [{Dq/ Dcrit) — 1] — Bgq. (8) 

For galaxies Dq 10 — 20. If we adopt Do/Dcrit ^ 2, km/kf ~ l/h, where I ~ lOOpc is 
the forcing scale of the turbulence, h ~ 400 — lOOOpc, then the mean field strenth would 
be 1/4 to 1/10 of equipartition at saturation. This estimate is more pessimistic than the 
more general limit Bt < {Q5y^'^{km/kfy^^Bcq, in section 2. This is basically because, 
for Dq exceeding but near Dcrit, the large-scale dynamo saturates even with modest a 
suppression, after which there is no further helicity transfer << 1). On the other hand, 
for Dq >> Dcrit Eq. (^) seems to suggest fields greater than obtained in section 2. This 
rather points to the limitations of the quasi-linear model for non-linear saturartion, in 
this case, than to a violation of the more general limit. 

A major caveat to the above limits is that galaxies have boundaries, and if one 
has a flux of helicity due to small-scale fields preferentially leaving the system then one 
may avoid the above constraints (cf. Blackman and Field 2000; Klccorin ct al. 2000). 
Artificial removal of small-scale fields in a simulation, periodically, does indeed lead to 
enhanced large-scale field growth; so the idea works in principle (Brandenburg, Dobler 
and Subramanian 2002). But so far the simulations which involve boundaries do not 
show a preferential out-flux of small-scale fleld helicity (Brandenburg and Dobler 2001). 
Another possibility is to think of a non- helical dynamo (Vishniac and Cho 2001), but 
there is no evidence yet for its working in a simulation by Arlt and Brandenburg (2001) 
designed to capture the effect. Clearly, thinking of ways out of the constraints implied 
by helicity conservation will be crucial to understand galactic magnetism. 
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